Macroscopic control parameter for avalanche models for bursty transport. 
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Similarity analysis is used to identify the control parameter Ra for the subset of avalanching 
systems that can exhibit Self- Organized Criticality (SOC). This parameter expresses the ratio of 
driving to dissipation. The transition to SOC, when the number of excited degrees of freedom 
is maximal, is found to occur when Ra — * 0. This is in the opposite sense to (Kolmogorov) 
turbulence, thus identifying a deep distinction between turbulence and SOC and suggesting an 
observable property that could distinguish them. A corollary of this similarity analysis is that SOC 
phenomenology, that is, power law scaling of avalanches, can persist for finite Ra, with the same 
Ra — * exponent, if the system supports a sufficiently large range of lengthscales; necessary for 
SOC to be a candidate for physical (Ra finite) systems. 



I. INTRODUCTION 

It is increasingly recognized that a large group of phys- 
ical systems can be characterized as driven, dissipative, 
out-of-equilibrium and having a conservation law or laws 
(see the comprehensive treatments of [HQ)- They usu- 
ally have many degrees of freedom (d.o.f.), or excited 
modes, and long range correlations leading to scaling or 
multiscaling. Two examples are fully developed turbu- 
lence (see e.g. 0, H|) and Self- Organized Criticality 

(SOC, 1,10). 

The SOC paradigm has found particular resonance 
with workers attempting to model, and understand, 
'bursty', scale free transport and energy release in mag- 
netized plasmas (for a recent review, see, for example 
8]). Simplified avalanche models have been proposed, 
and points of contact with the data investigated, in the 
astrophysical context; most notably to describe magne- 
tospheric activity ( d, [1(1 [HI G3 , see also the review |13j 
and refs. therein), the dynamics of the solar corona (e.g. 
[3, [ID, [TH, Q3| , see also the review [T^]), and accretion 
disks (e.g. [2(3, HH [Hj] ) . I n t ne context of magnetically 
confined laboratory plasmas there have been extensive 
efforts to construct avalanche models that make points 
of contact with the system under study, and to establish 
signatures characteristic of SOC dynamics in experiments 
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1371]). There have also been attempts to establish whether 
the signatures of SOC can emerge from magnetohydro- 
dynamic (MHD) or reduced MHD models, (e.g. in the 
solar coronal context [H, H§|). 

Since the original suggestion of Bak et al in 
that SOC "... could be considered a toy model of 
generalized turbulence" there has been continuing de- 
bate on the possible relationship of turbulence to SOC 
[4TI l42l l43l . 0, [4f| . Similarities in the statistical signa- 
tures of turbulence, and systems in SOC, have been noted 
(see e.g. [H|,[42]]). In particular, it has recently been ar- 
gued in the context of astrophysical plasmas that SOC 
and turbulence are aspects of a single underlying phys- 
ical process (see [48|, [49[ and references therein). How- 



ever, the extent to which SOC, as opposed to turbulence, 
uniquely captures the observed dynamics in magnetically 
confined laboratory plasmas (see [Hoi . [5ll l52j ) or is in- 
deed consistent with it (see [53, [54|) has been brought 
into question. Key observables such as power law distri- 
butions of patches are not unique to SOC systems (for 
an example see 1551 1 , also the comprehensive discussion in 

;n). 

Our focus here is then to establish the macroscopic 
similarities and differences between turbulence and SOC 
in the most general sense. A central idea in physics is 
that complex and otherwise intractable behavior may be 
quantified by a few measurable macroscopic control pa- 
rameters. In fluid turbulence, the Reynolds number Re 
expresses the ratio of driving to dissipation and param- 
eterizes the transition from laminar to turbulent flow. 
Control parameters such as the Reynolds number can 
be obtained from dimensional analysis (see e.g.0, [5r3|). 
without reference to the detailed dynamics. From this 
perspective the level of complexity resulting from the de- 
tailed dynamics is simply characterized by the number 
A of excited, coupled d.o.f. (or energy carrying modes) 
in the system. The transition from laminar to turbu- 
lent flow then corresponds to an (explosive) increase in 
A. The nature of this transition, the value of the Re at 
which it occurs, and the rate at which A grows with Re 
all depend on specific system phenomenology. Dimen- 
sional arguments, along with the assumptions of steady 
state and energy conservation, are however, sufficient to 
give the result that A" always grows with Re (as in [57l ]. 
see also 

We anticipate that an analogous control parameter for 
complexity, Ra, will exist for the wider group of systems 
discussed above. Interestingly it is now known that such 
a control parameter that expresses the ratio of driving 
to dissipation does indeed exist for SOC. In this paper 
we will give a prescription to obtain Ra generally from 
dimensional analysis, that is, without reference to the 
range of detailed and rich phenomenology that any given 
system will also exhibit. The rate at which A varies with 
Ra is again dependent on this detailed phenomenology. 
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We will see that similarity arguments, along with the 
assumptions of steady state and energy conservation, are 
however, sufficient to determine whether or not N grows 
with Ra- 

The question of control parameters in Self-Organized 
Criticality was initially controversial, as the name leads 
one to expect. It was originally argued (@, S3|, see also 
[58| ) that avalanching systems self organized to the SOC 
state without a tuning parameter. Subsequent analysis 
has established a consensus (see 0, 0, [H, MM, & that 
some tuning exists, at least in the sense that SOC is a 
limiting behavior in the driving rate h and the dissipation 
rate e, such that h/e — > with h,e — > 0, (and h < e, that 
is, a steady state). This understanding is exemplified in 
Jensen's constructive definition as given in [?J of SOC 
as the behavior of "slowly driven interaction dominated 
thresholded" (SDIDT) systems. Clearly then, h/e plays 
the role of a control parameter. This SDIDT limit h/e — » 
has been investigated extensively (e.g. [13, [HI), most 
recently with respect to finite size scaling in the limit of 
increasingly large system size (e.g. [62, 63]). 

Here we are concerned with the relevance to SOC to 
physically realizable systems, and in particular natural 
ones, where the system size is finite and the driving may 
be unknown and/or highly variable. Our focus is on pa- 
rameterizing the level of complexity of the system as we 
take it away from the SDIDT limit by increasing the 
driver, in a system of large but fixed size. For avalanche 
models exhibiting SOC, we will argue that distinct realiz- 
able avalanche sizes play the role of excited d.o.f. of the 
system. The SOC state is then characterized by max- 
imal excited d.o.f., that is, avalanches occurring on all 
lcngthscales supported by the system. Far from the SOC 
state, the system becomes ordered with few excited d.o.f. 
and exhibits laminar flow. The SDIDT limit is reached 
by taking Ra to zero and we will show that this indeed 
maximizes the number of excited d.o.f. N. The SDIDT 
limit is thus in the opposite sense to fluid turbulence 
which maximizes N at Re —> oo. 

This suggests a possible means to distinguish obser- 
vationally between turbulence and SOC in observations 
and experiments of driven, magnetically confined plas- 
mas. For example, a power law region in the power spec- 
tral density of some quantity that probes the flow is often 
identified in both laborator y a nd astrophysical confined 
magnetized plasmas (e.g. [13, [H, H^, [3(|) and is dis- 
cussed in the context of both SOC and turbulence. This 
power law region will always be of finite spatio- temporal 
range (an 'inertial range' of the cascade). Our results im- 
ply that this 'inertial range' will decrease as we increase 
the drive for SOC, whereas it will increase for turbulence- 
providing an experimental or observational test to distin- 
guish these phenomena. 

Our relationship between Ra and N implies the pos- 
sibility of large but finite N for small but non- zero Ra', 
hence an important corollary is that SOC phenomenol- 
ogy can quite generally persist under conditions of finite 
drive in a sufficiently large bandwidth system. This has 



been seen in specific avalanche models (see f&il l65l |66| ) . 
Here, since our result flows from dimensional analysis, 
we will see that this is a generic property of avalanching 
systems. 



II. SIMILARITY ANALYSIS AND CONTROL 
PARAMETER 

We shall focus on how the well-established techniques- 
similarity analysis, as described in 56|. and the H- 
theorem obtained by Buckingham in [67| nearly a cen- 
tury ago can be used to reveal new information about 
avalanche models exhibiting SOC (i.e. [ll.laldl7l.l58l.l59l]). 

The systems that we have in mind all have strongly 
coupled excited d.o.f. that transport some quantity from 
the driving to the dissipation scale. They have the fol- 
lowing properties: 

I. The many excited d.o.f. of the system are coupled; 
there is some dynamical quantity that freely flows over 
all the excited d.o.f. of the system. We can characterize 
a flux £i of this quantity associated with processes that 
occur on lengthscale I, that is, e; is the transfer rate of 
the dynamical quantity through I to neighboring length- 
scales. 

II. The system is not necessarily in equilibrium but is 
in a steady state on the average. 

III. The dynamical quantity is conserved so that given 
(II) the injection rate £j„j balances the dissipation rate 
Ediss, that is £i n j ~ Si ~ £diss in an ensemble averaged 
sense. 

IV. The solution is of a scaling type, that is: 

where a > and Lq and 51 are the largest and small- 
est lengthscales respectively that are supported by the 
system. 

V. The number of excited d.o.f. can be parameterized 
by a single macroscopic control parameter. 

We will identify the control parameter for these sys- 
tems in terms of known macroscopic variables by formal 
dimensional analysis (similarity analysis or Buckingham 
n theorem, see e.g. [5(1 [67j )■ The essential idea is that 
the system's behavior is captured by a general function 
F which only depends on the relevant variables Qi .v 
that describe the system. Since F is dimensionless it 
must be a function of the possible dimensionless group- 
ings, the Hi..m(Qi..v)i which can be formed from the 
Qi..v- The (unknown) function F(Jli, Tl2, ..ILm) is uni- 
versal, describing all systems that depend on the Qi..v 
through the TIi..m{Qi..v) an d the relationships between 
them. If one then has additional information about the 
system, such as a conservation property, the Hi..m(Qi..v) 
can be related to each other to make F explicit. Thus 
this method can lead to information about the solution 
of a class of systems where the governing equations are 
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TABLE I: II theorem applied to homogeneous turbulence. 
Variable dimension description 

Lo L driving length scale 

rj L dissipation length scale 

U LT^ 1 bulk (driving) flow speed 

v l?T~ x viscosity 



unavailable or intractable, often the case for complex sys- 
tems where there are a large number ( N here) of strongly 
coupled d.o.f.. If the V macroscopic variables are ex- 
pressed in W physical dimensions (i.e. mass, length, 
time) then there are M = V — W dimensionless group- 
ings. 

The properties I-V above restrict the choice of relevant 
Qi..v- First, we have only specified there is a transfer 
rate on lengthscale I, e\ (property I) of some dynamical 
quantity, its precise nature is irrelevant. Consequently, 
the only physical dimensions of the transfer rate £\ rele- 
vant to the problem are length and time, so that W = 2. 
Second, property (V) is that there is a single control pa- 
rameter, LTi which may be expressed as a function of 
the number of excited degrees of freedom N . To incor- 
porate the scaling property (IV) we will seek solutions 
such that n 2 = g(L /Sl) = f(N). This means that the 
system's behaviour is captured by some F(Hi,U.2) = C 
which fixes M = 2 (C is a constant). The LTi and II2 
are related to each other via properties II and III (con- 
servation and steady state). We then have that V — 4; 
there are always four relevant macroscopic variables to 
consider. 

T see this in action, we begin with a relatively well 
understood example, namely Kolmogorov (K41) turbu- 
lence. Our aim here is to straightforwardly illustrate the 
above approach by obtaining the control parameter, the 
Reynolds number Re as a function of N via dimensional 
analysis; for a detailed discussion of the universal scal- 
ing properties of K41 turbulence and their origin in the 
Navier Stokes equations see for example 3 . As above, for 
K41 we have four relevant macroscopic variables (given 
in Table 1) and two dimensionless groups: 

IL 1 = ^ = R E ,U 2 = ^ = f(N) (2) 

III is just the Reynolds number Re of the flow, and the 
ratio of lengthscales II2 is related to the number of d.o.f. 
N that can be excited. We now see how R E is related 
to f(N) by relating II ! to n 2 . For incompressible fluid 
turbulence, our dynamical quantity si is the time rate of 
energy transfer per unit mass through length scale I. The 
procedure is then as follows: 

1. conservation and steady state imply (ensemble av- 
eraged) that Si n j ~ £1 ~ Ediss', that is, the average 
energy injection rate £j n j balances the average en- 
ergy dissipation rate £di SS - 



2. the rate at which energy is transferred to the fluid 
is from dimensional analysis: Si n j ~ U 3 /L 

3. Dimensional analysis of Navier Stokes gives £diss ~ 

4. £inj ~ £<iiss then relates III to II2: 

(3) 

v \ 77 J 
and fixes exponent (3 = 4/3 

5. the solution is of scaling type, so that: 

(4) 

with a > by definition 

6. thus 

RE-^^-Nto (5) 

and (3 N = § > 

The value of the exponents a and /3 will depend on the de- 
tailed phenomenology of the turbulent flow. An estimate 
based on K41 for example, with (3 — 4/3 from the above, 
and a = D = 3 where D is Euclidean dimension Q, im- 
plies a high degree of disorganization and will be mod- 
ified for example if the turbulence is intermittent. Im- 
portantly, the only property of turbulence with which we 
are concerned here is that both (3 > and a > so that 
(3n = (3 1 a > 0. This identifies the Reynolds number as 
the control parameter for a process (turbulence) which 
simply excites more active modes or d.o.f. as we increase 
Re- 

We now see how above arguments apply to other sys- 
tems as defined above, in particular to avalanche models. 
Without recourse to details of the system, similarity anal- 
ysis will be sufficient to obtain the relationship between 
the control parameter R and the number of degrees of 
freedom N of the form: 

R - N@ N (6) 

The value of the exponent (3n will depend on the details 
of these systems but crucially we will see that the sign of 
Pn is fixed by the similarity analysis. This is sufficient 
to establish whether or not, as in the case of turbulence, 
increasing R increases the the number of excited d.o.f. in 
the system. 

III. CONTROL PARAMETER FOR 
AVALANCHING SYSTEMS 

We now envisage a generic avalanche model in a sys- 
tem of size Lq where the height of sand is specified on 
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TABLE II: II theorem applied to an avalanching system. The 
sand carries a property with dimension S. 



Variable dimension description 


L L 


system size 




51 L 


grid size 




e ST' 1 


system average 


dissipation/loss rate 


h ST' 1 


average driving 


rate per node 



a grid, with nodes at spacing 51. Sand is added to in- 
dividual nodes, that is, on length scale 51 at an average 
time rate Ei n j = h per node. There is some process, here 
avalanches, which then transports this dynamical quan- 
tity (the sand) though structures on intermediate length 
scales SI < I < La- Sand is then lost to the system (dis- 
sipated) at a time rate e over the system size Lq. On 
intermediate length scales 51 < I < Lq, sand is conserva- 
tively transported via avalanches (see also [H, [69|, l70l]). 
In our discussion here we follow [5( and assume that the 
transport timescale is fast, so that avalanches occur in- 
stantaneously and do not overlap. There must be some 
detail of the internal evolution of the pile that maximizes 
the number of length scales I on which avalanches can oc- 
cur. For avalanche models this is the property that trans- 
port can only occur locally if some local critical gradient 
is exceeded; as a consequence the pile evolves through 
many metastable states. If these length scales represent 
excited d.o.f. then the number N of d.o.f. available will 
be bounded by Lq and SI so that N ~ (Lo/5l) a , with 
D > a > for D > 1 (a may be fractional). 

The four relevant variables for the avalanching system 
are given in Table 2. The two dimensionless groups are: 



III 



Ra, n 2 



Lo 
51 



f(N) 



(7) 



We will now relate the control parameter LTi = h/e 
to the number of excited d.o.f. by following the same 
procedure as above, £/ now refers to the time rate of 
transfer of 'sand' through length scale I. 

1. conservation and steady state imply (ensemble av- 
eraged) e irl j ~ El ~ £ d iss 

2. In Euclidean dimension D there are (Lq/SI) d 
nodes; D > by definition. The rate at which 
'sand' is transferred to the pile is then from dimen- 
sional analysis: £i n j ~ h{L$/5l) D 

3. the system average dissipation rate is defined as 

£ — E-diss 

4. Ei n j ~ £diss then gives h(L /5l) D ~ e or: 



R 



h 



51 



D 



< \L 

thus in the above notation fixes (3 — — D < 



(8) 



5. the number TV of d.o.f. available will be bounded 
by L and 51 so that: 



AT ~ -if 



(9) 



with D > a > for D > 1 (the value of a depends 
on the details and may be fractional). 



6. thus 



and j3 



N 



h 




e 




(± = . 

a 


-s. <o 

a 



N- 



(10) 



We then have that the number of excited d.o.f. decreases 
as we increase the control parameter Ra — h/e. Thus 
we recover the SDIDT limit for SOC, namely Ra — > 0, 
but now explicitly identify this limit with maximizing 
the number of excited d.o.f.. Our result from dimen- 
sional analysis is to obtain Ra ~ N@ N and to show quite 
generally that for the avalanching system (3n < 0. 

Our dimensional analysis for the avalanche model maps 
onto that for K41 turbulence, so in that sense Ra = Re, 
that is, Ra is the avalanching system's 'effective Reynolds 
number', which expresses the ratio of driving to dissi- 
pation. Both Re and Ra increase with driving of the 
system, but the system's response is quite different. In 
the case of K41 turbulence, the system can excite more 
modes or degrees of freedom and the flow becomes more 
disorganized, whereas in the avalanche models, less d.o.f. 
are available so the system is pushed toward order. The 
essential difference between the two systems in this con- 
text is as follows. As we increase the driving in K41 
turbulence, the smallest lengthscale rj can decrease (via 
Navier Stokes) to provide the necessary dissipation to 
maintain a steady state, and since we have assumed scal- 
ing the system simply excites more modes or d.o.f.. On 
the other hand, in the avalanche models both the smallest 
and largest lengthscales are fixed; increasing the driving 
will ultimately introduce sand at a rate that exceeds the 
rate at which sand can be transported by the smallest 
avalanches, as we discuss next. 



IV. SOC- LIKE BEHAVIOUR UNDER 
INTERMEDIATE DRIVE 

For avalanching to be the dominant mode of trans- 
port of sand, there are conditions on the microscopic de- 
tails of the system; specifically, there must be a separa- 
tion of timescales such that the relaxation time for the 
avalanches must be short compared to the time taken 
for the driving to accumulate sufficient sand locally to 
trigger an avalanche. Avalanches are triggered when a 
critical value for the local gradient is exceeded, the criti- 
cal gradient can be a random variable but provided it has 
a defined average value g, we have that on average, we 
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would need to add gSl sand to a single cell of an initially 
flat pile to trigger redistribution of sand. The number of 
timesteps that this would take to occur would on aver- 
age be (gSl) I (hSt) where again 51 is the cell size and ft is 
the timestep. This gives the condition for avalanching to 
dominate transport on all lengthscales in the grid [51, Lo], 
so that avalanches only occur after many grains of sand 
have been added to any given cell in the pile and is the 
strict SDIDT^iH limit: 

hSt < g6l (11) 

We will now consider an intermediate behaviour (see also 
0) 

g5l < hSt < gSl (jj^j (12) 

where the driver is large enough to swamp of order 
h6t/(gSl) cells of the pile at each timestep (each addition 
of sand), but this is still much smaller than the largest 
avalanches that the system is able to support since the 
largest possible avalanche in a system of Euclidean di- 
mension D is (Lo/6l) D cells. 

For a given physical realization of the sandpile, that 
is, fixed box size L and grid size 61, successively in- 
creasing h5t above g6l then successively increases the 
smallest avalanche size (to some 6V > 51). Ultimately 
as h and hence Ra is increased to the point where 
hSt ~ gSl (Lo/6l) D there will be a crossover to laminar 
flow, as each addition of sand drives avalanches that are 
on the size of the system. 

We now assume that the avalanching process is self- 
similar, so that the system is large enough that the prob- 
ability density of avalanche sizes 5 is P(S) ~ 5 -7 over a 
large range of 5; that is, finite sized effects do not domi- 
nate. Consequently this intermediate, finite Ra behavior 
will be 'SOC like', with avalanches occurring within the 
range of lengthscales [51', Lq] with power law statistics 
sharing the same exponent 7 as at the SDIDT limit. 

We will illustrate these remarks with simulations of the 
Bak, Tang and Wiesenfeld (BTW) (see [H) sandpile in 
2D, where the driving occurs randomly in time and is 
spatially restricted to the 'top' of the pile. In all cases 
shown, the critical gradient (threshold for avalanching) 
is g = 45l~ 1 , and normalized distributions of the number 
of topplings in an avalanche 5 are shown (we take top- 
plings as a measure of avalanche size following [f|). In 
Figure 1 we plot the results from two simulations in a 
box of size Lq/51 = 100, under driving rates h = 4<5t — 1 
and h! = 16ft" 1 . We can see that as we increase the 
driving rate from h = 4ft -1 to h' — 16<5£ — 1 , the occur- 
rence probability of the smallest avalanches is reduced, 
and on these normalized histograms, the probability of 
larger events is increased. These larger events for the run 
with ti = 16ft" 1 , that is, for 5 ~ [10 2 - 10 3 5 ], still fol- 
low the same power law scaling as the h = 4ft -1 run (the 
precise location on the plot of the crossover in behaviour 
will depend on details of the dynamics of the pile) . This 



is to be anticipated provided that transport on these in- 
termediate scales is still dominated by avalanching, that 
is, intermediate scale avalanches still have the property 
that they relax on a timescale that is much faster than 
that required by the driving to initiate an avalanche. If 
this is the case, then the phenomenology of these inter- 
mediate scale avalanches is unchanged by the increase in 
the driving rate and as a consequence, except close to 
the crossover in statistics, their scaling exponent is, as 
we see, unchanged. As the system has self similar spatial 
scaling we can also anticipate obtaining the same solu- 
tion for these avalanches subject to a rescaling; 5, which 
is a measure of avalanche size, will simply scale with hSt, 
the sand which must be redistributed at each timestep 
since h5t > g5l . This is shown in the lower plot of Figure 
1 where we have rescaled the h! = 16ft -1 intermediate 
range driving results by 5 — > 5/16. We can see that 
power law regions of the plots that both correspond to 
avalanching now coincide. 

We can go further and anticipate that two realizations 
of the system, one with h and Lq and the other with h! = 
Ah and L' Q = ALq give the same solution for P(S) under 
rescaling 5 — > S' /A D . This is shown in Figure 2 where 
we compare two runs of the sandpile (i) with h = 4ft _1 
and L = 10051 and (ii) with h' = 16 and L' = 40061, 
in the same format as Figure 1. We can indeed see a 
close correspondence of the avalanche statistics in the 
power law region of the plot once we have rescaled the 
b! = 16ft -1 and L' = 40051 run by 5 -> 5/16 (at the 
largest 5, the histograms do not precisely collapse under 
this self affine scaling, see [10] for a discussion of the finite 
size type of scaling properties of the model) . 

This establishes a general property of avalanching sys- 
tems that has been seen in several representative SOC 
models, such as in [H, [65|, [6(| . Depending on the details, 
specifically, provided that a separation of timescales for 
avalanching can be maintained, some SOC systems will 
show scaling in systems where the drive is in fact highly 
variable. One could argue that such robustness against 
fluctuations in the driving is necessary for SOC to pro- 
vide a 'working model' in real physical systems where the 
idealized SDIDT limit may not be realized. 



V. CONCLUSIONS 

We have used similarity and dimensional analysis to 
discuss high dimensional, driven, dissipating, out of equi- 
librium systems, in particular avalanching systems that 
exhibit bursty transport that can be in an SOC state. 
These act to transport a dynamical quantity (e.g. for 
the avalanche models, sand) from the driving to the dis- 
sipation scale, in a manner that is conservative, that is 
steady state on the average, and that shows scaling. The 
generic nature of this method of analysis implies that 
our results are not restricted to sandpile models per se, 
and have wider application to physical systems that show 
bursty transport, and scaling. We have postulated that a 
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FIG. 1: (Color online) Avalanche size normalized distribu- 
tions for two runs of the 2D BTW[1, [U sandpile driven at 
the top corner formed by two adjacent closed boundaries, the 
other boundaries are open. Lo/Sl — 100 and h5t = 4 (•) and 
hSt — 16 (x); (left) probability densities; (right) as (a) with 
probability density for the h8t = 16 avalanche sizes rescaled 
S -» 5/16. 



'class' of these systems have a single control parameter R 
which expresses the ratio of the driving to the dissipation 
and which can be related to the number of excited de- 
grees of freedom N. Dimensional analysis then leads to a 
relationship of the form R ~ N@ N and without reference 
to any detailed phenomenology of the system, determines 
the sign of (3n- 

We have focussed on avalanche models that can exhibit 
SOC, for which the above identifies the control parame- 
ter Ra = h/e. The limit Ra — > is just the well known 
SDIDT limit of SOC. Specific avalanching systems will 
have different values of (3n but will all share the essen- 
tial property that we obtain here, that /3jv < so that 
that N is maximal under the limit of vanishing driv- 
ing. Our formalism for SOC has close correspondence 
with that for Kolmogorov homogeneous isotropic turbu- 
lence. A minimalist interpretation of our results is that 
Kolmogorov turbulence maximizes the number of excited 
d.o.f. N under maximal (infinite) driving in contrast to 
SOC. A maximalist interpretation is that Ra is analogous 
to the Reynolds number Re- This establishes an essen- 
tial distinction between turbulence and SOC. Practically 
speaking, it can for example arise because if we fix the 



10 10 10 10 10" 



10 

s 



FIG. 2: (Color online) Avalanche size normalized distribu- 
tions for Lo/Sl = 100, HSt = 4 (•) and ^ = 400, hSt = 16 
(+); (left) probability densities; (right) as (a) with probability 
density for the hSt = 16 avalanche sizes rescaled S — » .5/16. 



outer, driving scale in Kolmogorov turbulence, the dissi- 
pation scale can simply adjust as we increase the driving. 
Since the system shows scaling, this acts to increase the 
available degrees of freedom. Avalanching on the other 
hand, is realized in a finite sized domain (box) and driven 
on a fixed, smallest scale, so increasing the driving be- 
yond a certain point simply swamps the smallest spatial 
scales, thus reducing the available degrees of freedom. In- 
creasing the driving then pushes Kolmogorov turbulence 
toward increasingly disorganized flow, and avalanching 
systems toward more ordered (laminar) flow. 

A corollary is that SOC phenomenology, that is, power 
law scaling of avalanches, can persist for finite Ra with 
the same exponent that is seen at the Ra — > limit, 
provided the system supports a sufficiently large range 
of lengthscales. This has been seen previously for spe- 
cific realizations of avalanche models [64| but is shown 
here to be quite generic; and is a necessary property for 
SOC to be a candidate for physical (Ra finite) systems. 
As the driving is increased, the excited number of de- 
grees of freedom (modes) decreases for SOC and increases 
for turbulence, so that in principle one could distinguish 
SOC from turbulence observationally by testing how the 
bandwidth (range of spatio- temporal scales) over which 
scaling is observed, varies with the driving rate. 



7 



Acknowledgments the EPSRC, and the NSF (under grant No. NSF PHY05- 

51164). 

We thank M. P. Freeman and K. Rypdal for discus- 
sions. This research was supported in part by the STFC, 



[1] D. Sornette, Critical Phenomena in Natural Sciences, 2nd 
erf., Springer, Berlin, (2004). 

[2] J. P. Sethna, Statistical Mechanics- Entropy, Order Pa- 
rameters and Complexity, Oxford University Press, Ox- 
ford, (2006). 

[3] U. Frisch Turbulence. The legacy of A.N. Kolmogorov, 
Cambridge University Press, Cambridge, (1995). 

[4] T. Bohr, M. H. Jensen, G. Paladin, A. Vulpiani, Dynami- 
cal Systems Approach to Turbulence: Cambridge Nonlin- 
ear Science Series (No. 8), Cambridge University Press, 
Cambridge, (1998). 

[5] P. Bak, C. Tang and K. Wiesenfeld Phys. Rev. Lett 59, 
381 (1987). 

[6] V. Frette, K. Christensen, A. Malthe-Sorenssen, J. Feder, 

T. Jossang and P. Meakin Nature, 379, 49 (1996). 
[7] H. Jensen Self organised criticality: Emergent complex 
behaviour in physical and biological systems, Cambridge 
University Press, Cambridge, (1998). 
[8] R. O. Dendy, S. C. Chapman, M. Paczuski, Plasma Phys. 

Cont. Fusion, 49, A95, (2007) 
[9] T. S. Chang, IEEE Trans. Plasma Sci. 20, 691 (1992). 
[10] Chapman S. C, Watkins N. W., Dendy R. O., Helander 
P. and Rowlands G., Geophys. Res. Lett. 25, 2397, (1998) 
[11] Lui A. T. Y., Chapman S. C, Liou K., Newell P. T., 
Meng C.-L, Brittnacher M. J. and Parks G. K., Geophys. 
Res. Lett., 27, 911, (2000) 
[12] Uritsky V. M., Klimas A. J. Vassiliadis D., Geophys. Res. 

Lett, 30, 1813 doi:10.1029/2002GL016556 (2003) 
[13] S. C. Chapman, N. W. Watkins, Space Science Re- 
views special issue "Challenge to Long-Standing Un- 
solved Space Physics Problems in the 20th Century", 
95,293, (2001). 

[14] Lu E. T. and Hamilton R. J. Astrophys. J., L89, 380, 
(1991) 

[15] D. W. Longcope, E. J. Noonan, Astrophys. J., 542, 1088, 
(2000) 

[16] A. Bershadskii, K. R. Sreenivasan, Eur. Phys. J. B, 35, 
513, (2003) 

[17] D. Hughes, M. Paczuski, R. O. Dendy, P. Helander, and 
K. G. McClements, Phys. Rev. Lett, 90, 131101 (2003) 

[18] L. Morales, P. Charbonneau, Astrophys. J., 682, 654, 
(2008) 

[19] P. Charbonneau, S. W. Mcintosh, H.-L Liu, T. J. Bog- 
dan, Solar Phys., 203, 321, (2001) 

[20] S. Mineshige, M. Takeuchi, and H. Nishimori, Astrophys. 
J., 435, L125 (1994). 

[21] K. M. Leighly and P. T. OBrien, Astrophys. J., 481, L15 
(1997). 

[22] R. O. Dendy, P. Helander, M. Tagger, Astron. Astro- 
phys., 337, 962 (1998). 

[23] P. H. Diamond, T. S. Hahm, Phys. Plasmas, 2, 3640 
(1995) 

[24] Newman D. E., Carreras B. A., Diamond P. H. Hahm T. 

S. (1996) Phys. Plasmas 3 1858, (1996) 
[25] B. A. Carreras, D. Newman, V. E. Lynch, P. H. Diamond, 



Phys. Plasmas, 3, 2903 (1996). 
[26] B. A. Carreras et al., Phys. Rev. Lett, 80, 4438 (1998). 
[27] X. Garbet and R. Waltz, Phys. Plasmas, 5, 2836 (1998). 
[28] Y. Sarazin and P. Ghendrih, Phys. Plasmas, 5, 4214 

(1998). 

[29] T. L. Rhodes, R. A. Moyer, R. Groebner, E. J. Doyle, R. 
Lehmer, W. A. Peebles, C. L. Rettig, Phys. Lett. A, 253, 
181 (1999). 

[30] P. A. Politzer, Phys. Rev. Lett, 84, 1192 (2000). 

[31] Chapman S. C, Dendy R. O. and Hnat B., Phys. Rev. 

Lett. 86 2814 (2001) 
[32] S. C. Chapman, R.O. Dendy, B. Hnat, Phys. Plasmas, 8, 

1969, (2001) 

[33] R. Sanchez, B. Ph. van Milligen, D. E. Newman, B. A. 

Carreras, Phys. Rev. Lett, 90, 185005 (2003) 
[34] T. K. March, S. C. Chapman, R.O. Dendy, J. A. Merri- 

fleld, Phys. Plasmas, 11, 659, (2004) 
[35] S. C. Chapman, R. O. Dendy, B. Hnat, Plasma Phys. 

Cont. Fusion, 45,301, (2003) 
[36] L. Garcia, B. A. Carreras, D. E. Newman, Phys. Plasmas, 

9, 841 (2002) 

[37] L. Garcia, B. A. Carreras, Phys. Plasmas, 12, 092305, 
(2005) 

[38] L. Vlahos, M. K. Georgoulis, Astrophys. J., 603, L61, 
(2004) 

[39] G. Einaudi, M. Velli, H. Politano, A. Pouquet, Astrophys. 

J., 457, L113, (1996) 
[40] P. Bak, C. Tang and K. Wiesenfeld Phys. Rev. A 38, 364 

(1988). 

[41] M. Paczuski and P. Bak Phys. Rev. E, 48, R3214 (1993). 
[42] G. Boffetta, V. Carbone, P. Giuliani, P. Veltri, and A. 

Vulpiani Phys. Rev. Lett. 83, 4662 (1999). 
[43] S. T. Bramwell et al. Phys. Rev. Lett, 84, 3744 (2000) 
[44] N.W. Watkins, S.C. Chapman, and G. Rowlands, Phys. 

Rev. Lett, 89, 208901, (2002). 
[45] S. C. Chapman, G. Rowlands, N. W. Watkins, J. Phys. 

A. Math. Gen. 38, 2289, (2005). 
[46] M. De Menech, A. L. Stella, Physica A, 309, 289, (2002) 
[47] K. R. Sreeenivasan, A. Bershadskii, J. J. Niemela, Phys- 
ica A, 340, 574, (2004) 
[48] M. Paczuski, S. Boettcher, M. Baiesi, Phys. Rev. Lett, 

95, 181102, (2005) 
[49] V. M. Uritsky, M. Paczuski, J. M. Davila, S. I. Jones, 

Phys. Rev. Lett, 99, 025001,(2007) 
[50] J. A. Krommes, Phys. Plasmas, 7, 1752, (2000) 
[51] J. A. Krommes and M. Ottaviani, Phys. Plasmas, 6, 3731 

(1999) 

[52] P. A. Politzer, M. E. Austin, M. Gilmore, G. R. McKee, 
T. L. Rhodes, C. X. Yu, E. J. Doyle, T. E. Evans, R. A. 
Moyere, Phys. Plasmas 9, 1962 (2002) 

[53] E. Spada, V. Carbone, R. Cavazzana, L. Fattorini, G. 
Regnoli, N. Vianello, V. Antoni, E. Martines, G. Seri- 
anni, M. Spolaore, and L. Tramontin, Phys. Rev. Lett, 
86, 3032 (2001) 

[54] V. Antoni, V. Carbone, R. Cavazzanal, G. Regnoli, N. 



8 



Vianellol, E. Spada, L. Fattorini,2, E. Martines, G. Se- 

rianni, M. Spolaore, L. Tramontin, P. Veltri, Phys. Rev. 

Lett, 87, 045001 (2001) 
[55] S C Chapman, R O Dendy and N W Watkins, Plasma 

Phys. Cont. Fusion, 46, B157, (2004) 
[56] G. I. Barenblatt Scaling, self- similarity, and intermediate 

asymptotics, Cambridge University Press, Cambridge, 

(1996). 

[57] A. N. Kolmogorov, C. R. Acad. Sci., 30, 301, (1941) 
[58] J. P. Sethna, K. A. Dahmen and C. R. Myers Nature, 

410, 242 (2001). 
[59] R. Dickman, M. A. Munoz, A. Vespignani, and S. Zapperi 

Brazilian J. Phys., 30, 27 (2000). 
[60] A. Vespignani and S. Zapperi Phys. Rev. E, 57, 6345 

(1998). 

[61] M. Vergeles, A. Maritan, J. R. Banavar Phys. Rev. E, 

55, 1998(1997). 
[62] G. Pruessner, O. Peters, Phys. Rev. E, 73, 025106(R), 

(2006) 



[63] M. J. Alava, L. Laurson, A. Vespignani, S. Zapperi, Phys. 

Rev. E, 77, 048101 (2008), G. Pruessner, O. Peters, Phys. 

Rev. E, 77, 048102 (2008) 
[64] N. W. Watkins, S. C. Chapman, R. O. Dendy, P. 

Helander, G. Rowlands Geophys. Res. Lett., 26, 2617 

(1999). 

[65] A. Corral and M. Paczuski Phys. Rev. Lett, 83, 572 
(1999). 

[66] V. M. Uritsky, A. J. Klimas and D. Vassiliadis Phys. Rev. 

E, 65, 046113 (2002). 
[67] Buckingham, E., Phys Rev., 4, 345, (1914). 
[68] J. X. de Carvalho and C. P. C. Prado Phys. Rev. Lett, 

84, 4006 (2000), 
[69] K. Christensen, D. Hamon, H. J. Jensen, and S. Lise 

Phys. Rev. Lett, 87, 039801 (2001) 
[70] J. X. de Carvalho and C. P. C. Prado, Phys. Rev. Lett, 

87, 039802 (2001). 



